Order of the chiral phase transtion in the a models 
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Abstract 

The chiral phase transition at a certain critical temperature is a restoration mechanism of the chiral symmetry, 
broken by the nonzero mass of quarks and mesons. The transition can be studied through several models, among 
which are the a models. An analysis is made on the linear and nonlinear a models with different approximations, 
and we show that the transition is of second order in both cases. 



1 Introduction 

A chiral transformation is a rotation in the isospin space 
of quarks. When a quark system is invariant under such 
a transformation, we say that the system is chirally sym- 
metric. This is true in the case of massless quarks, but we 
know that actually they have mass. So with regard to the 
light quarks only (up and down), we can say that the chiral 
symmetry is an approximated symmetry of the strong in- 
teractions jH]. Nonzero mass of the corresponding mesons 
also breaks the symmetry. 

There are two kinds of chiral transformations. The par- 
tially conserved Noether current is a vector in one case, 
and an axial vector in the other case. When working with 
the up and down quarks, the transformations can be stud- 
ied through the group SU(2)xSU(2), which is isomorphic 
to 0(4). This isomorphism is the reason for the Heisen- 
berg magnet model being the unique description of the low 
energy dynamics of the QCD. 

The a models arise in this context. The a model has 
the usual kinetic term and the potential 
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with a field 3> having N components, where A is a positive 
coupling constant, and / w is the pion decay constant; this 
model is renormalizable. In the limit that A — ► oo the po- 
tential goes over to a ^-function constraint on the length 
of the field vector. This is the so-called nonlinear a model. 
Both models describe a quark system in which the chiral 
symmetry is broken, in consistency with the nonzero mass 
of the light quarks and/or the light mesons. 

If quarks are massless, QCD is expected to undergo a 
chiral phase transition, which may have implications for 
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high energy nucleus- nucleus collisions. A pion gas with a a 
meson undergoes the chiral phase transition. In this work, 
we use both the a models to describe it and derive its pres- 
sure. The linear a model is studied through a mean field 
approximation, while the nonlinear a model is studied us- 
ing chiral perturbation theory The first and second term of 
the perturbation series are considered. In both studies, the 
transition is of second-order and the critical temperature is 
the same, as we will show later. 

The developments towards the calculation of the criti- 
cal temperature follow closely those of Bochkarev and Ka- 
pustaPj. 



2 Linear sigma model 



The linear a model Lagrangian is 
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where A is a positive coupling constant. The bosonic field 
<I> has N components. Rather arbitrarily, we define that 
they are N — 1 ir fields and a a field as the iVth com- 
ponent. The space spanned by these components has the 
usual Euclidean metric. Since the a has the quantum num- 
ber of the vacuum, and since the symmetry is broken at low 
temperatures, we immediately allow for a a condensate v 
whose value depends on the temperature 1/(3 and is yet to 
be determined. 
We write 



$i(x,t) = 7Tj(x, t), i=l,.. 

$jv(x,i) = v + a(x,t). 



,N-l; 



(3) 
(4) 
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In terms of these fields the Lagrangian becomes 

% = i {d^f + I (d,af — j (v 2 — flf - | (« 2 - f 2 ) (2va + a 2 + ir 2 ) — j (2w + a 2 + n 2 ) 2 



(5) 



We define the effective masses 



Reorganizing terms, 

2 = -\ (jl v 2 ) 2 + \ \{d^f ~ rhln 2 + (d,af - m 2 ^ 2 ! - Xv (v 2 - f 2 ) a - Xva (a 2 + n 2 ) - \ (a 2 + n 2 ) 2 (6) 



At zero temperature, the potential is minimized when 
v = f„. In this particular case, 



ml = X {v 2 - f 2 ) = 0, 
ml = \(3v 2 -f 2 )=2\f 2 , 



(7) 

(8) 



so the Goldstone theorem is satisfied. 

The action at finite temperature is obtained by rotat- 
ing to imaginary time, r = it, and integrating r from to 
j3 = l/T. We keep the Minkowsi metric: = d/dx* 1 with 
do = d/dt = id I dr. Then the action is 



S = -±(f 2 -v 2 ) 2 pV 



+ {d,a) 2 - m 2 a a 2 - Xv (v 2 - f 2 ) a 



— Xva [a 2 
where the abbreviation 
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dr / d 3 a: 
Jv 
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is used. 

At any temperature v is choosen such that (a) = 0. 
This eliminates any one-particle reducible diagrams in per- 
turbation theory, leaving only one-particle irreducible (1PI) 
diagrams. We will allow for v to be temperature depen- 
dent. The simplest approximation at finite temperature 



is the mean field approximation; in this, we will neglect 
interactions among the particles or collective excitations. 

The pressure of a free relativistic boson gas can be writ- 
ten in several ways: 



P (T,m) 



ln(l 



(11) 

(12) 



i r d 3 k 

d 3 fc k 2 1_ 

(2tt) 3 3uj - 1 

The pressure of the pion gas is given by 

P = -jy^m = -±(f 2 -v 2 f + VPo(T 7 m a ) 

+ V(N-l)P (T,m w ) (13) 

The mean field approximation thus calculated has been 
a very simple, very good approximation, allowing to ana- 
lyze lots of phenomena in finite temperature field theory. 
It was used in all the pioneering papers. 

2.1 second-order phase transition 

We expect that the temperature rise tends to disorder the 
condensate until it dissapears. In the linear a model we 
analyze this phenomena by expanding the free boson gas 
pressure about zero mass: 



P (T,m) 
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The efective masses are proportional to the square root of A, and terms with m 3 or A 3 / 2 are not part of the mean 
field approximation, so we neglect the last term of the expansion (|14fl . The pion gas pressure becomes 



P{T,v)^-\{fl~v 2 ) 2 
X 



k 2 4 T 2 2 \ , x /tt 2 , T 2 „ 
90 T * "24"* J+^W "24^ 



N ^ - o - 2 f 2 -^T 2 
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(15) 
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Now we maximize the pressure with respect to v: 
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N+2 T 2 
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The pressure must be a maximum because we expect 
the pion gas to be in a stable or mctastable state. We 



verify the value of the second derivative for both solutions. 



d 2 P 
dv 2 



N ■ 
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with respect to the cases of Eq. Ijl8l) : the first case repre- 
sents a sure maximum. 

The condensate goes to zero at a critical temperature 
given by 



T 
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Above this temperature the condensate is zero, because 
thermal fluctuations are very large. 
We evaluate Eq. l|2"0)l again: 



d 2 P 
dv 2 



,2 _ N+2 



N + 2 (rp2 _ rp2 

12 \ c J 



A i^(T 2 -T 2 ) 



v 2 = 0. 



(22) 



In the second case the pressure is a minimum at T < T c , 
so we drop it here. 



Now we will show that this is a second-order phase transition, calculating the pressure under the two temperature 
regimes. 



P<(T) =N— T 4 + -V 
90 2 
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/ 2 _ T 2y + '2fl ?L 2 _ N_ 2 I _ ^ 2 2 



(-2 fOT^ _ rp2\ _ ^ f2rp2 

48 -M^ ± c) 24 /7r ' 
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(23) 
(24) 
(25) 
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Evaluating the pressure in the critical temperature, 



PAT) = N—T 
<v ' 90 ' 



48 ^ c ~ 24 J " c ~ 90 c 



XN 
48" 



/ 2 T 2 = P>(T) 



(28) 
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Evaluating the first derivative, 
P' <{ T) = N^T 3 + A (^±2) 2 (T 2 - T 2 ) T 
P^T)=N 2 -^T 3 + ^LflT 
P'ATo) = N 2 -^T? + ^f 2 T c = P^{T C ). 

Evaluating the second derivative. 

P»(T) = N^-T 2 + A (^±2) 2 (3T 2 - T 2 ) 



XN 
~12 



(29) 
(30) 
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,2tt 2 
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P»(T)=N^T 2 + — ft 

P i'(T c ) = N^T 2 + X {^f 2T 2 + ^Lf 2 
= P ;(T C ) + X(^) 2 2T 2 ; 
we find that it is discontinuous. 

3 Nonlinear sigma model 

This model is defined by the Lagrangian 

^ = i(^*) 2 , 
together with the constraint 

/* = <£ 2 0M) 



(31) 



(32) 



(33) 



(34) 



The length of the chiral field is fixed and cannot be changed 
by thermal fluctuations, so we say that chiral symmetry 
is built into this model and, therefore, there can be no 
chiral-symmctry restoring phase transition. Taking the 
limit A — > oo constraints the length of the chiral field to 
be / w just as in the nonlinear model. The critical tempera- 
ture, however, is independent of A at least in the mean field 
approximation. So it would seem that the phase transition 
survives. If it is true, then we ought to be able to derive 
it entirely within the context of the nonlinear a model. 
That is what we shall do, although it involves a lot more 
effort than the treatment of the linear model. Since the 
only parameter in the model is f^, and we are interested in 
temperatures comparable to it, we cannot do an expansion 
in powers of T/ f w . The only parameter is N, the number 
of field components. This suggests an expansion in powers 
of 1/JV. 

We begin by representing the field-constraining 5 func- 
tion by an integral. 



m@b' cxp / d 4 z (jSf + ib' ($ 2 - f 2 )) . (35) 



The zero frequency, zero momentum condensate v is al- 
lowed. Following Polyakov, we also separate the zero fre- 
quency, zero momentum mode of the auxiliary field b' 



Integrating over all the other modes will give us an ef- 
fective action involving the constant part of the fields. We 
will then minimize the free energy with respect to these 
constant parts, which is a saddle point approximation. In- 
tegrating over fluctuations about the saddle point is a finite 
volume correction and of no consequence in the thermody- 
namic limit. The Fourier expansions are 

$*(x,t) =7r i (x,r) = ^^e l ( x -P+^^(p,n), (36) 

p , 71 

$jv(x,t) =w + ct(x,t) (37) 

= v + \[v EP'" el(x ' p+ "" r) ^P' n )- ( 38 ) 

(39) 



6'(x,r) = i— + &(x,t) 



V ^ 



b{p,n), 



where we define the abbreviation 



(40) 



(41) 



The zero frequency and zero momentum modes have 
been excluded from the summations. The field $ must 
be periodic in imaginary time for the usual reasons but 
there is no such requirement on 6, hence Lu n — 2imT and 
v n = ■nnT . Since b has dimensions of inverse length squared 
we inserted another factor of T so as to make its Fourier 
amplitude dimensionless, as it is for the other fields. From 
(|35p. the action becomes 



S = 



d 4 z(^ [{d^) 2 + (d,{ V + a)f} 

+ l (^ + b)^ 2 + {v + a) 2 -f 2 )^ (42) 
J d 4 x U [{d^f - mV + (d^a) 2 - m V] 
- ib(2va + 7T 2 + a 2 )) + \ra 2 (fl - v 2 )(3V. (43) 



Please note that terms linear in the fields have been inte- 
grated to zero because 



(tt 4 ) = (a) = (b) = 0. 



(44) 



An effective action is derived by expanding e s in pow- 
ers of b. The term linear in b vanishes n account of 
6(0,0) oc (b) = 0. The term prooprtional to b 2 is not zero 
and is exponentiated, thus summing a whole series of con- 
tributions. The term proportional to b 3 is not zero either 
and may also be exponentiated, summing an infinite series 
of higher-order terms left out of the order b 2 exponenti- 
ation. After the expansion is made, we perform a path 
integration over the pion and a fields, on ly on the term 
quadratic in b 2 . At last, the terms are scaled: b — > b/\2N. 
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We have 



S E = - i Yl + P 2 + m2 ) [*KPj ») ' *(-P. ~ n ) 

II(p, oj n ,T,m) 



p.n 



+ a(p,n)a(-p,-n) --^ Yl 
2 v 2 



+ 



N to 2 + p 2 + m 2 



b(p, 2n)6(-p, -2n) 



±m 2 (/ 2 - i; 2 )/3V + ^(& 3 /ViV). 



(45) 



The effective action is an infitine series in 6. The co- 
efficients are frequency and momentum dependent, arising 
from one-loop diagrams. In addition, each successive term 



is supressed by 1/yN compared to the previous one. This 
is the large N expansion. 

The propagators for the pion and a fields are of the 
usual form 



D/(p,w„,m) = w 2 +p + m 2 



and the propagator for b is 



D b 1 (p,w n ,m) = U(p,uj n ,T,m) + 



2 v 2 

N U) 2 + p 2 + TO 2 



(46) 



(47) 



with m and v to be determined. There appears the one- 
loop function 



n( P ,c n ,r,m) = r^|^(K-^) 2 + (p-fc) 2 + m 2 ) 1 (w? + A 2 + m a ) *) 



(48) 



Keeping only the terms up to order b 2 in Se (The rest vanishes when — > 00) allows us to obtain an explicit expression 
for the partition function and the pressure. This includes the next-to-leading order in N: 



Z = e W{fl-jyv exp f_ NVT 2^ J- In (/3 2 ( 
xexp Ll^T 2 ^ln 



u>l+p 2 + m 2 )) 



U(p,uj n ,T,m) 



N u! 2 + p 2 + in 2 



P=^nZ= \m 2 (f 2 v 2 ) - £ T £ In (/3 2 ( W 2 + P 2 + m 2 )) | ^ In 



6(p,2n)6(-p,-2n) 
2 



n(p,w„,T, m) 



iV w 2 + p 2 + m 2 



(49) 
(50) 



The second term under the last logarithm should and will be set to zero at this order. It may be needed at higher 
order in the 1/N theory to regulate infrared divergences. 
The pressure is extremized with respect to m 2 : 



dP T 



dm 2 VZ 
_T 1 
~V2 



J 2><$>®b' exp [J d 4 x [jz? + ib' ($ 2 - f 2 )]^j i J d 4 a- [-tt 

(-7T 2 - ( T 2 + / 2 - U 2 ) = -(7r 2 )-( ( T 2 )+/ 2 - U 2 =0 



2 a 2 + f 2 -v 2 ]=0 



(51) 
(52) 



This condition is equivalent to the thermal average of 
the constraint: 



f 2 =v 2 + (n 2 ) + (a 2 ) = (* 2 ) 



(53) 



If an approximation to the exact partition function is 
made, such as the large TV expansion, this constraint should 
still be satisfied. In fact, it will isolate a preferred value for 
m. 

To leading order in N we may neglect the term involving 
II entirely. The pressure is then 



P = ^m 2 (f 2 -v 2 )+NP (T,m). 

This must be a maximum with respect to v: 

DP 2 

— — = — rn v = 0, 

ov 



(54) 



(55) 



which is equivalent to the condition (a) = by means of 
l)4()[l . There are two possibilities: 



1. m = 0: There are massless particles, or Goldstonc 
bosons, and the value of the condensate v is deter- 
mined by the thermally averaged constraint. This is 
the asymmetric phase. 



2. v = 0: The thermally averaged constraint 153|) is 
satisfied by a nonzero T-dependent mass. There 
are no Goldstone bosons. This is the symmetric, or 
symmetry-restored phase. 



Evidently, there is a chiral symmetry-restoring phase tran- 
sition. 
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In the leading order of the 1/N approximation, the par- 
ticles are represented by free fields with a potentially In- 
dependent mass m. For any free bosonic field <j>, 



dP 
dm 2 



d 3 p 1 1 



(2tt) 3 wePv-l' 



(56) 



with lo 2 = p 2 + m 2 . Thus, extremizing the pressure with 
respect to m 2 is equivalent to satisfying the thermally av- 
eraged constraint: 







dm 2 



p,7Z 



9 



1 



1 T V 1 

~~ ~9 L^i ,,,2 i 



2 e —' + p 2 + m 

p,n " 

2 (2tt) 3 2w V e ' - 1 



We renormalize the expression: 



2~ ^CJ 2 +W 2 
p,n ,L 



(57) 



1 T / d 3 p P 1 



(2tt) 3 



,3" - 1 



Then we have 



dm 2 



(T, m) 



<* 2 ) 



(58) 



(59) 



•••U4^-^-y^ 2 > = (60) 

,. £ =, 2 + <. 2 ) + (a 2 ) = (* 2 ) (!) (61) 
We consider now the two different phases: 

• The mass is zero in the asymmetric phase. The con- 
straint is satisfied by a temperature-dependent con- 
densate: 

v 2 {T) = fl-^ (62) 

This condensate goes to zero at a critical temperature 
of 

1 c ~ jfU- 



(63) 



At exactly Tc, the thermally averaged constraint is 
satisfied by the fluctuations of N massless degrees of 
freedom without the help of a condensate. 

• In the symmetric phase the condensate is zero. The 
constraint is satisfied by thermal fluctuations alone: 



fl = N(&)=N 



d 3 p 1 1 



(2?r) 3 



(64) 



At fixed T, thermal fluctuations decrease with increas- 
ing mass. The constraint is only satisfied by massless exci- 
tations at one single temperature, namely Tc- At temper- 
atures T > Tc the mass must be greater than zero. Near 



Tc the mass should be small, and the fluctuations may be 
expanded about m = as 



fl ~NT 
12 



1 

12 



m 



4ttT 



&ir 2 T 2 
rn 



In 



4ttT 16tt 2 T 2 



(65) 



3mT 



(66) 



As T approaches Tc from above, the mass approaches 
zero like 

m = h ^ T ^ ■ (67) 

This is a second-order phase transition since there is no 
possibility of metastable supercooled or superheated states. 
We will evaluate the order of this transition in the cases 



a) m 

b) v = 



0. 



0. 



Function evaluation Given m > 1 in the neighborhood 
of Tc we have for any T, and in particular, for T = Tc, 



P a (T) =NP (T, 0) = A^T 4 



(68) 



T 2 f f 2 + NP (T,^(T 2 



-TP 



= (T - Tc 

24 9 v c 



T 



T 2 1 n J 



V90 



24 2 9T 2 " 1 ~ c 



=N*-T* + ^-(T 
90 24 9 v 



T c ) 



T 



1 N ir 2 , , 
(T 2 

2 24 9 v 



First derivative evaluation 

r 2 



n'(T< 



Ar 90 4T ^ = ^ (Tc) 



Second derivative evaluation 

Nn 2 



=p:(t 



4 {ZTl-T^ 

24 9 V c c 



Att 2 

y"9 



(69) 
(70) 

(71) 

(72) 

(73) 
(74) 



This is a second order phase transition. 

The mass must grow faster than the temperature at 
very high temperatures in order to keep the field fluctua- 
tions fixed and equal to f 2 . Asymptotically, the particles 
move nonrelativistically. This allows us to compute the 
fluctuations analytically. We get 



A 



3/2 



me 



-m/T 



(75) 
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This is a trascendental equation for m. IT can also be 
written as 



=Tln 



NT mT 



2ttU V ^ft 



T 2 

Tin -j (76) 
1 c 



It is rather amusing that, in leading order of the large 
N approximation, the elementary excitations are massless 
below Tq, become massive above Tc, and at asymptotically 
high temperatures move nonrelativistically. 

The result to first order in l/N provides a good insight 
into the nature of the two-phase structure of the nonlinear 
a model, but it is not quite satisfactory for two reasons: 

1. It predicts N massless Goldstone Bosons in the bro- 
ken symmetry phase when we know there ought to be 
2V-1. 



2. In this model, T c 



12f 2 /N, and in the linear a 
model it is T c = 12f 2 /(N + 2). We expect them to 
be the same in the limit A — > oo. 

Both these problems can be rectified by inclusion of the 
second order in l/N; this is, the inclusion of the b field. 

It is natural to expect that the b field will contribute 
essentially one negative degree of freedom to the T 4 term 
in the pressure so as to give TV — 1 Gioldstonc bosons in 
the low temperature phase. Therefore, we move one of the 
N degrees of freedom and put it together with the b con- 
tribution. 



= l-MZ=\m 2 (fl-v 2 ) 



V 2 

N - 1 



T 



p.n 



p,7l 

oil, ,2,2, 2\ 



Fiti+r +m 2 )U(p,iu ni T,m) 



(77) 



After the procedure described by Bochkarev and Ka- 
pusta^Q, we do an expansion in m/T as before. 



P = (N 1)^T 4 - ^m 2 T 2 + lm 2 (f 2 - v 2 ) + £_m*T 

(78) 

We maximize in the high temperature phase where 
v = 0: 



W±2 mT 2 +m/ 2 + A 

f 2 =^T 2 - £ m T 

_rp2 \ N+2 _ N_m] 
1 I 12 4ir T J ■ 



(80) 
(81) 



This gives the same critical temperature as in the mean 
field treatment of the linear a model: 



Tr 



12 



N + 2 



J 7T 



(82) 



The mass approaches zero from above like 

m(T) = ^ ( f Ar t 2) (T 2 - T c) (83) 



We will evaluate the order of this transition in this order 
of the xPT, in the cases 

a) m = 0, 

b) v = 0. 

Function evaluation Given m > 1 in the neighborhood 
of Tc we have for any T, and in particular, for T = Tc, 



P a (T) =(N - 1)P (T, 0) = (N — 1)— T 4 

P b (T) =(N - 1)— T 4 - — (jV+ ? 3 (T 2 - T 2 
v ' 1 y 90 24 9iV 2 v 



(84) 



c) 



TT 2 (N + 2) 2 



(T 2 - Tc) f 2 



18N 2 T 2 
n 2 {N + 2f 2 2 3 
27 12T 2 N 2 K c) 

First derivative evaluation 

Pi(T c )^(N -l)^-AT c = P' a {T c ) 
Second derivative evaluation 



Pl'{T c ) =N—12T t 



7T 2 (N + 2f 



4 3r-T, 



c 



(85) 



(86) 



(87) 



90' c 24 9iV 2 

n 2 (N + 2) 3 2 
P a (Tc)- Y -^^T c 



This is a second order phase transition. 

4 Concluding remarks 

Actually, the chiral phase transition has been studied us- 
ing three approximations: Mean field approximation on the 
linear a model, first order of xPT, and first and second or- 
der of xPT. We have shown in the three approximations 
that the transition is of second order, However, the critical 
temperatures are not the same, as pointed out for the first 
order of xPT. This made the inclusion of the second order 
necessary. 
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